Astronomy & Astrophysics manuscript no. aa 


©ESO 2013 


March 1,2013 





(N 

o 

(N 



O 
< 

o 

6 

o3 



> 

O 
O 
O 

o 

(N 



13 



The magnetic field of IRAS 16293-2422 as traced by shock-induced 

H 2 masers 

F. 0. Alves 1 , W. H. T. Vlemmings 2 , J. M. Girart 3 and J. M. Toixelles 4 

1 Argelander-Institut fur Astronomie, University of Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
e-mail: falves@astro.uni-bonn.de 

2 Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden 
e-mail: wouter . vlemmings@chalmers . se 

3 Institut de Ciencies de l'Espai (IEEC-CSIC), Campus UAB, Facultat de Ciencies, C5 par 2 a , 08193 Bellaterra, Catalunya, Spain 
e-mail: girart@ice . cat 

4 Institut de Ciencies de l'Espai (CSIC)-UB/IEEC, Universitat de Barcelona, Marti i Franques 1, E-08028 Barcelona, Spain 
e-mail: torrelles@ieec . cat 

Received ; accepted 

ABSTRACT 

Context. Shock-induced fLO masers are important magnetic field tracers at very high density gas. Water masers are found in both 
high- and low-mass star-forming regions, acting as a powerful tool to compare magnetic field morphologies in both mass regimes. 
Aims. In this paper, we show one of the first magnetic field determinations in the low-mass protostellar core IRAS 16293-2422 at 
volume densities as high as 10 8 ~ 10 cnT 3 . Our goal is to discern if the collapsing regime of this source is controlled by magnetic fields 
or other factors like turbulence. 

Methods. We used the Very Large Array (VLA) to carry out spectro-polarimetric observations in the 22 GHz Zeeman emission of 
H 2 masers. From the Stokes V line profile, we can estimate the magnetic field strength in the dense regions around the protostar. 
Results. A blend of at least three maser features can be inferred from our relatively high spatial resolution data set (~ 0.1"), which 
is reproduced in a clear non-Gaussian line profile. The emission is very stable in polarization fraction and position angle across the 
channels. The maser spots are aligned with some components of the complex outflow configuration of IRAS 16293-2422, and they 
are excited in zones of compressed gas produced by shocks. The post-shock particle density is in the range of 1 - 3 x 10 9 cnT 3 , 
consistent with typical water masers pumping densities. Zeeman emission is produced by a very strong line-of-sight magnetic field 
(£~113mG) 

Conclusions. The magnetic field pressure derived from our data is comparable to the ram pressure of the outflow dynamics. This 
indicates that the magnetic field is energetically important in the dynamical evolution of IRAS 16293-2422. 

Key words, stars: formation - masers - polarization - ISM: magnetic fields - ISM: individual objects: IRAS 16293-2422 

1 . Introduction suit regarding how magnetic fields can be resolved in very small 

spatial scales was obtained by Vlemmings et al. (2006), who car- 
Spectro-polarimetric observations of masers are a powerful ried QUt Very Long Baseline ^ (VLBA) polarization obser- 
tool to study the magnetic field properties in the maser pump- yations of H Q maserfj arQund the c heus A HW2 hi h . mass 
ing zone. Water masers are unique because they are found in a object 
variety of astrophysical environments and, in particular, in star- 
forming regions at distinct mass regimes. In contrast, methanol Thermal dust polarization emission is difficult to observe in 
and OH masers are found associated mostly with high-mass star the most embedded portions of molecular clouds due to the low 
formation sites (e.g., Vlemmings 2008; Sanna et al. 2010a,b). sensitivity to polarized data (and the feasibility is then limited 
The most commonly observed water maser line is the (6 16 - 5^) toward a handful list of objects). Submillimeter (submm) emis- 
transition at 22 GHz, an excellent probe of molecular gas at very sion suffers from depolarization effects like unresolved fields or 
high volume densities (n H , in the 10 s to 10 10 cnT 3 range, Elitzur roundness of dust grains (Goodman et al. 1992, 1995; Lazanan 
et al. 1989). While the Zeeman splitting is small (~ 10 3 Hz et al - 1997 )- Water masers can be thought as an observational 
yuG 1 ), the line strength is sometimes sufficient for measuring tGo1 to overcome this issue, since they are excited at very high 
the circular polarization and therefore for a direct measurement densities. 

of the line-of-sight (LOS) component of the magnetic field. If IRAS 16293-2422 (hereafter, 116293), is a prototypical low- 
linear polarization is also measured, then the full 3-D magnetic mass protostellar system located in the p Ophiuchus molecu- 
field configuration can be derived. lar cloud (d =* 120 pc, Loinard et al. 2008). This source is a 
Recently, several studies have appeared in the literature re- well-studied binary system, usually referred as source A and B, 
porting the efficiency of maser spectro-polarimetry technique us- with angular separation of 5" (600 AU) and very embedded in a 
ing masers. These studies have revealed the magnetic field prop- dense molecular core (Wootten 1989; Looney et al. 2000). Both 
erties in very dense molecular environment around the circum- sources have a very rich chemistry which is typically found in 
stellar envelopes of evolved stars and young stars (e.g., Surcis hot cores (Ceccarelli et al. 2000; Kuan et al. 2004; Bisschop et al. 
et al. 201 la,b; Perez-Sanchez et al. 201 1). A representative re- 2008; J0rgensen et al. 201 1). Although component B appears to 
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be a single source, high resolution interferometer data reveal a 
higher degree of fragmentation for source A. In fact, two dust 
submm components were detected within this object with an an- 
gular separation of ~ 0.6" (~ 72 AU, Chandler et al. 2005). 
116293 has two large scale (~ 1') bipolar CO outflows, one of 
them associated with source A while the powering source of the 
other CO outflow is a matter of debate (Walker et al. 1988; Stark 
et al. 2004; Yeh et al. 2008). Recently, observations of the SiO 
(8-7) emission have revealed a compact molecular outflow also 
associated with source A (Rao et al. 2009). 116293 has strong 
water maser emission that has been well monitored (Wilking & 
Claussen 1987; Terebey et al. 1992; Claussen et al. 1996; Furuya 
et al. 2003). The strongest features appear at blueshifted and red- 
shifted LSR velocities with respect to the ambient cloud velocity 
(~ 4 km s~'), typically between Vlsr of -5 and 10 km s _1 . The 
H2O maser emission very often has intensities of more than 100 
Jy, and at times it is stronger than 300 Jy, at LSR velocities of 
~ 7 km s -1 . 

Tamura et al. (1993) performed observations of the 1.1 mm 
dust polarized emission toward 116293 at an angular resolution 
of 19". They found that the magnetic field lines are perpendicu- 
lar to the major axis of the dense elongated disk-like molecular 
structure (~ 10" size, ~ 1200 AU) observed in C ls O by Mundy 
et al. (1990). This configuration was corroborated recently by 
Rao et al. (2009), who used the Submillimeter Array (SMA) and 
obtained a polarization map at an angular resolution of ~ 2" 
(~ 240 AU), resolving both sources A and B. The mean volume 
density traced by the SMA maps is ~ 6x 10 7 cirT 3 . The polariza- 
tion pattern around source A is compatible with a hourglass mor- 
phology for the magnetic field, whose strength was estimated in 
~ 4.5 mG. The SMA maps show that there is a misalignment be- 
tween the outflow direction and the magnetic field axis, roughly 
in agreement with model predictions where the magnetic energy 
is comparable to the centrifugal energy. In contrast, source B is 
associated with an uniform and apparently undisturbed magnetic 
field. 

In this work, we report spectro-polarimetric Very Large 
Array (VLA) observations of FLO masers toward 116293. In § 2, 
we describe the details of the observational setup. In § 3, the re- 
sults obtained from the maser spectroscopy are shown. In § 4, 
we discuss a possible correlation between the masers features 
and the outflows. The magnetic field strength estimation and its 
implication on the core dynamical evolution are discussed in sec- 
tions 5 and 6, respectively. Finally in § 7 we summarize our con- 
clusions. 

2. Observations 

The observations were done with the VLA (NRAO 1 , New 
Mexico, USA) in its extended A-configuration, on 2007 June 
25 th and 27 th . The tracks lasted ~ 5.5 hours each. A total of 27 
antennas were used, 10 of them already retrofitted with the new 
system, resulting in a combined VLA/EVLA (Extended VLA) 
observation. We used the K band receivers (22-24 GHz) tuned 
at the frequency of the water maser (6 16 - 5i3) rotational tran- 
sition (vo = 22.23508 GHz). We used the full polarization ca- 
pability of the correlator, selecting a bandwidth of 0.7813 MHz 
(~ 10.5 km s in velocity). The spectral setup contains a total 
of 128 channels which covers most of the velocity range of the 
strongest water maser features (observed around the brightest 

1 The National Radio Astronomy Observatory - NRAO - is a facility 
of the National Science Foundation, operated under cooperative agree- 
ment by Associated Universities, Inc. 



feature previously reported at ~ 7 km s~') at the spectral resolu- 
tion of 0.08 km s~ 1 . The quasar J 1 626-298 was used as gain cal- 
ibrator. The quasar J1331+305 and the radio source J1751+096 
were used for polarization calibrations in order to obtain cor- 
rections in the instrumental feed polarization and the corrected 
position angle of the polarization vectors. All three calibrators 
were also used for bandpass corrections. We performed self- 
calibration using the channel with the strongest intensity. Since 
the detected emission is unresolved, both phase and amplitude 
self-calibration solutions were applied to the other channels. 

The combined VLA/EVLA setup introduces about 8% clo- 
sure errors on VLA-EVLA baselines-based corrections on the 
bandpass. This happens due to the distinct bandpass response of 
the EVLA antennas compared to the VLA ones. These errors 
become more severe at narrow bandwidths, when the process- 
ing of the digital signal of EVLA antennas aliases the power 
response at the band edges. During the data reduction, we took 
special care on determining the bandpass solutions by perform- 
ing particular calibration strategies recommended to the user in 
the NRAO webpage 2 . They consisted of flagging EVLA-EVLA 
baselines prior to bandpass calibration and unflagging them prior 
to applying the bandpass solutions to the data set. The bandpass 
corrections are then applied to an averaged single-channel multi- 
source dataset ("channel 0" file) and to the multi-channel origi- 
nal file. 

Data reduction was done with the Astronomical Image 
Processing Software package (AIPS). Imaging of Stokes param- 
eters I, Q, and U were generated with a quasi-uniform weight- 
ing (robust = -1, Briggs 1995). Maps of polarized fraction (P) 
and position of polarization angles (PA) were obtained by com- 
bining Stokes Q and U images in such a way that P — y = 

and PA = i tan '(^). The resulting synthesized beam is 

0. 14" x 0.08", with a position angle of -5.7°. The rms noise for 
channels where no emission is detected is ~ 8 mJy beam 1 , and 
it increases up to 23 mJy beam -1 at the peak emission channel. 
A slightly lower rms is observed for polarized intensity. 

3. Results 

3.1. H 2 maser line 

The contour channel maps of the H2O emission observed 
with the VLA/EVLA data are shown in Fig. 1. The emission ex- 
tends over a wide range of velocities (4.5 < Vlsr < 9 km s _1 ), 

1. e., at redshifted LSR velocities with respect to the cloud sys- 
temic velocity (~ 4 km s _1 ). Our spectral setup covers only a 
small portion of LSR velocities blueward of it (~ 2 < Vlsr < 4 
km s~') but no emission was detected within this spectral range. 
The water maser line has peak intensity of 170 Jy beam -1 , de- 
tected at Vlsr - 7.4 km s _1 . This emission was extensively re- 
ported in previous surveys also as the brightest feature (see, for 
instance, Claussen et al. 1996; Furuya et al. 2003). Some chan- 
nels have emission extended slightly eastward (Vlsr between 

6.5 and 7 km s _1 ), suggesting that the observed emission in 
these channels is marginally resolved. 

A scheme of the distribution of dust and molecular outflows 
in 116293 is shown in Fig. 2. The masers detected in our obser- 
vations are associated with source A. They are located ~ 0.25" 
(~ 30 AU in projection) to the southeast of the dust conden- 
sation Aa resolved by the subarcsecond submillimeter (submm) 
observations of Chandler et al. (2005). 

2 http : //www . via . nrao . edu/astro/guides/ evlareturn/ 
postproc/index . shtml#closure- line 
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Fig. 1. Channel maps of the H 2 emission toward IRAS 16293-2422. Contours are -50, 50, 500, 3xl0 3 , lxlO 4 , 2xl0 4 times 8 mjy beam" 1 , the 
rms noise of the maps with lower intensities (see Section 2). The peak flux of 167 Jy beam -1 occurs at V LSR = 7.4 km s~'. The LSR velocity of 
each channel is labelled. 



The Stokes I spectrum of the maser emission is shown in 
Fig. 3 (upper panel). The non-Gaussian line profile indicates that 
there are unresolved components in the spectrum. Apart of the 
peak intensity at 7.4 km s~', unresolved emission seems to be 
also present at lower velocities (Vlsr — 5.5 km s~' ) with a strong 
flux of ~ 20 Jy beam -1 . Fainter emission (~ 5 Jy beam -1 ) is 
also observed at higher velocity channels (Vlsr - 9.2 km s -1 ). 
Therefore, there are at least three unresolved components (Table 
1). A two-dimensional Gaussian profile fit on each of those 
components provides a mean spatial separation of ~ 22 milli- 
arcseconds (mas), which is higher than the precision on the rela- 
tive position determination (rms ; , OJ 2 mas, estimated from the 
ratio between the width of the synthesized beam and the signal- 
to-noise ratio of the fainter component). 

3.2. Polarized emission 

The spectrum of linearly polarized intensity is very similar 
to the Stokes I line profile except for the flux scale. It peaks at 
the same systemic velocity as the Stokes I spectrum, although 
it is weaker by a factor of ~ 30. The dependence of the linear 
polarization intensity, polarization fraction, position angle and 
Stokes I with systematic velocity is shown in Fig. 3. The mea- 
sured polarization fraction is 2.5 ± 0.2%. The polarization po- 
sition angle 6 is -23° and shows only small changes across the 
maser {erg = 2°), implying that the polarization vectors at differ- 
ent velocities trace basically the same region. The formal uncer- 
tainty in PA {erg = jjf-^f-, Wardle & Kronberg 1974) is small 
and ranges between ~ 0. 14° and 0.80° considering channels with 



Table 1. H 2 maser components in IRAS 16293-2422 



Vlsr 
(kms -1 ) 


^peak 

(Jy beam - ') 


a (2000) 
(h m s) 


6 (2000) 
(o / //) 


5.7 


23 


16 32 22.8830 


-24 28 36.495 


7.4 


168 


16 32 22.8808 


-24 28 36.487 


9.2 


5 


16 32 22.8818 


-24 28 36.493 



a Equatorial coordinates derived with the JMFIT task of AIPS. 



strong and weak linear polarization emission, respectively. The 
linear polarization as derived from Stokes U and Q maps is ex- 
hibited in Fig. 4, which shows the distribution of polarization 
vectors in the brightest velocity channel. The peak of polarized 
intensity is offset 0.05" with respect to the Stokes I peak. The 
polarization vectors can be parallel or perpendicular to magnetic 
field orientation in the plane of the sky. In § 5, we discuss which 
assumptions must be considered in order to solve this ambiguity. 

The line profile of the circular polarization (Stokes V) has 
the characteristic S-shape (see Fig. 5, lower panel). The Stokes 
V spectrum is proportional to the first derivative of the Stokes I 
spectrum and to the LOS component of the magnetic field. The 
red dashed line in Fig. 5 indicates the fit which best represents 
this derivation. Due to remaining gain differences among the two 
polarizations, scaled down replicas of Stokes I spectrum might 
contaminate the Stokes V data (Sarma et al. 2001; Vlemmings 
et al. 2002). These features were considered in the fit and sub- 
tracted from the synthetic V spectrum. The fraction of circular 
polarization, calculated as (V max - V m j n )/I max , is ~ 0.45% for the 
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Fig. 2. Scheme of the distribution of dust and molecular material in 
116293. The plus signal indicates the position of the peak intensity of 
our VLA water maser data. The ellipse is the deconvolved size of the 
dust continuum source A as derived by Rao et al. (2009). Triangles 
denote the position of the submillimeter condensations observed by 
Chandler et al. (2005). Stars denote the VLBI water maser detections 
of Imai et al. (2007) and straight lines denotes the direction of the CO 
and SiO outflows associated with this core (Rao et al. 2009). 



brightest component. The remaining maser features show only 
residual Zeeman profiles with amplitudes at the rms level. 

4. H 2 masers and the SiO/CO outflows 

Fig. 6 shows that the three possible maser features unre- 
solved in our data are distributed linearly. The velocity shift 
observed between Vlsr —5.7 km s _1 and higher velocities is 
oriented in a E-W direction (PA 110°). The proximity with 
source Aa and the complex outflow configuration observed in 
this zone suggest that the maser emission is being pumped in the 
dense circumstellar material around this object. 

Rao et al. (2009) report a strong and very young SiO outflow 
extending toward the northwest-southeast (NW-SE) direction, 
with the SE lobe detected at redshifted LSR velocities (see Fig. 
2). This SE component is also traced by a CO (3-2) outflow at 
V LSR ^ 9 km s 1 (Fig. 7 in Rao et al. 2009). The CO outflow 
has a dominant E-W orientation and has already been reported 
in the literature (Yeh et al. 2008). Rao et al. (2009) claim that 
the SiO and CO outflows are centered on sources Ab and Aa, 
respectively. Nevertheless, it is hard to disentangle the powering 
sources from their molecules maps, since that the separation be- 
tween Aa and Ab (~ 0.6") is much smaller than the HPBW of 
their SiO and CO maps (~ 3"). 

Our maser features are also detected at redshifted velocities 
and lie in a similar direction as the CO outflow. On the other 
hand, their positions south of Aa suggest that they may be ex- 
cited by the SiO outflow. In fact, the SE lobes observed for both 
SiO and CO at redshifted velocities (Fig. 6 and 7 from Rao et al. 
2009) match quite well to the LSR of our highest velocity fea- 
tures (Vlsr ~ 9 km s~'). In addition, VLA 3.7 cm data from 
Chandler et al. (2005) show two sources, Al and A2, close to 
Aa (see Fig. 6). These authors refer to A2 as a protostar which 
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Fig. 3. Stokes I (black line) and linear polarization intensity (multiplied 
by a factor of 10, red dashed line) spectra of the water maser emis- 
sion (upper panel). The spectra of polarization fraction (black line, left 
scale) and position angle (dashed line, right scale) are also shown (lower 
panel). 



powers a radio jet due to the bipolar shape seen in their high- 
resolution maps. On the other hand, source Al would be an ion- 
ized region produced by shocks between this flow and nearby 
dense gas. Therefore, it is possible that the water masers are gen- 
erated in the interaction regions traced by the radio source s and 
the outflows. 

Water masers detections were also reported toward source 
A through very long baseline (VLBI) observations (Imai et al. 
2007). The milliarcsecond resolution data revealed a spot exactly 
over Aa and another one to the SW of it (see Fig. 2). While the 
former may have been excited in the circumstellar gas of source 
Aa, the latter is likely created by shocks in the E-W outflow. 
The redshifted lobe of this emission has an open shell structure, 
therefore masers off from the outflow main axis are expected. 

5. The line-of-sight magnetic field strength in 116293 

From the Zeeman splitting formalism, the magnetic field 
strength can be correlated to the fraction of circular polarization 
by (see, for instance, Fiebig & Guesten 1989) 



= 2 x A p-pi x (Bcos6»)/AV/, 



(1) 

(2) 



where the A f _f , coefficient depends on the maser rotational lev- 
els F and F' , the intrinsic thermal linewidth Av,;, and the maser 
saturation degree, while Avj is the FWHM of the total power 
spectrum. For the Af-p coefficient, we adopted the range of 
0.012 to 0.018, since these values are consistent with models 
and observations (Nedoluha & Watson 1992; Vlemmings et al. 
2002). The water maser spectrum is a blend of different velocity 
components, so we adopt a linewidth ranging from 0.75 km s _1 , 
a typical value found in other regions (Vlemmings et al. 2006), 
to 1.0 km s -1 , the maximum estimated value in 116293 from our 



4 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



-24 28 36.1 



36.2 



z 
2 
t- 



36.3 



36.4 



36.5 



36.6 



36.7 



36.8 



36.9 



16 32 22.91 




22.90 



22.89 22.88 22.87 22.86 
RIGHT ASCENSION 2000) 



22.85 



Fig. 4. Distribution of H2O linear polarization vectors in the brightest 
emission channel (V LSR = 7.4 km s _1 ). Contours are -50, -30, 30, 50, 
500, 3 x 10\ 1 x 10 4 , 2 x 10 4 x 8 mjy beam 1 . Only polarization 



vectors whose P 
beam. 



1% are plotted. Vectors are sampled as 1/2 of a 



spectral profile. For these values, we find that Bios ranges from 

94 to 188 mG. The negative signal is inferred from the 

Stokes V shape, meaning that the field is pointing toward the ob- 
server. For narrower linewidths (~ 0.75 km s _1 ), which is more 
realistic for resolved water maser lines, the field strength ranges 

from 94 to 141 mG (see Table 2). We note, however, that 

given that the linear polarization fraction is about 3%, this maser 
is likely unsaturated and the field strength determination is a fair 
approximation. 

Fig. 4 shows the linear polarization vectors at the velocity 
channel with the brightest emission. There is a degeneracy be- 
tween the position angle of the linear polarization vectors and 
the magnetic field direction projected in the plane-of-sky (POS). 
The position angle of the polarization will be parallel or per- 
pendicular to the magnetic field in the POS for 9 > or 
9 < 9 cri , = 55°, respectively (Goldreich et al. 1973), where 9 
is the angle between B and the maser propagation direction and 
9 cr it is the so called Van Vleck angle. Given that the linear polar- 
ization fraction P/ depends on 9 and the maser saturation level, 
we could in principle constrain the value of 9 from our data and 
solve for this ambiguity. We used a non-LTE radiative transfer 
model in order to fit the observed intensity / and polarization 
fraction P/. However, since our line is a blend of features spa- 
tially unresolved, we were unable to determine 9. Although ob- 
servations at more extended configurations are needed to achieve 
this goal, we can at least claim that the POS field topology is 
quite ordered in both cases, i. e., it would remain ordered if ro- 
tated by 90°. 

According to the submm data of Rao et al. (2009), the di- 
rection of the dust polarization vectors associated with source A 
averages to 95°, which implies that the POS magnetic field lines 
have a PA of 5° with respect to the North direction. Assuming 
that the 22 GHz polarization vectors in Fig. 4 are parallel to the 
POS magnetic field, the main direction of the field lines is -23° 
(or its supplement, 157°) with respect to the North direction, as 
reported in Section 3.2. On the other hand, if our polarization 
vectors are perpendicular to the POS field, their main direction 
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Fig. 5. Stokes I (upper panel) and Stokes V (lower panel) spectra of the 
water maser emission. The red dashed line is the scaled derivative of 
total power / in order to guide the eye. 

would then be 67°. Although two distinct scales are being com- 
pared, it is interesting to point out that in both cases there is 
clearly a change in the direction of the field lines in the small 
scales traced by our VLA data. In the case of parallel field lines 
and polarization vectors, such a change could be produced by 
the interaction of the SiO outflow and the magnetic field, given 
that the field lines lie almost in the same direction of the outflow 
(PA ~ 145°, Fig. 2). 

6. Preshock magnetic fields and post-shock 
densities 

The only previous estimation of the magnetic field strength 
toward 116293 was performed by Claussen et al. (2003), who 

used VLBI observations of water masers to estimate Bios 40 

mG. Their maser spots appear associated with the Al cm source 
(Fig. 6), with a projected separation of ~ 20 AU at Vlsr ~ 2.5 
km s _1 (Claussen, private communication). Since this LSR ve- 
locity is similar to the core ambient velocity, it is unlike that this 
emission traces the same region observed with our VLA data. 
In any case, our work is then one of the first determinations of 
the field strength in a low-mass young stellar object at densities 
larger than 10 8 cm 4 . As discussed in § 4, the observed maser 
emission is excited in zones of compressed gas produced by 
shocks between the outflows and the ambient gas. In an ionized 
medium, shocks compress the gas by a factor tha, the Alfvenic 
Mach number defined as vJva, where v s is the shock velocity 
and \>a the Alfven speed. The magnetic field Bo in the pre-shock 
zone is also enhanced by a factor niA in the shocked gas so that 
B s ~ niA x Bo- Since va varies with Bo(n) -1 ^ 2 , then B s is propor- 
tional to v s (n) ' 2 , where n is the volume density of the pre-shock 
region. Rearranging this expression, the pre-shock density can 
be written as (Kaufman & Neufeld 1996): 



no 



ImG/ lkms-W 



cm 
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Fig. 6. Possible independent maser features (stars) as derived by gaus- 
sian fits. The accuracy in the position of each maser feature is ~ 2 mas, 
while the positions of the submm and cm sources are accurate within 
~ 25 mas and ~ 10 mas, respectively. The numeric labels are the maser 
velocities with respect to the Local Standard of Rest (LSR). The 3.7 cm 
sources Al and A2 (triangles) and the submm sources Aa and Ab (cir- 
cles) are also shown (Chandler et al. 2005). The sizes of each symbol 
are not proportional to their deconvolved sizes or their intensities. 

Table 2. S-field strength calculation from Eq. 2 (assuming AV/ = 0.75 
km s -1 ) 

A F - F > = 0.012 A F . F , = 0.015 A F . F , = 0.018 - 



Slos (mG) 



■141 



-113 



-94 



The LOS magnetic field strength calculated in § 5 can be 
rescaled in order to obtain the total field strength B lo , = B s 
in the shocked zone. According to Crutcher et al. (1999), the 
squared LOS magnetic field (Blos) 2 averaged over several lines 
of sight in the Galaxy is ~ B 2 ot /3. The shock velocities are re- 
trieved from submm water emission detected in shocked zones 
around 116293 (Ristorcelli et al. 2005). These data were modeled 
as C-shocks (low velocity, non-dissociative shocks, Kaufman & 
Neufeld 1996) with shock velocities between 12-15 km s 4 . 
Applying those values to equation 3, we obtain pre-shock den- 
sities of ~ 4.2*2 4 x 10 8 cm -3 , which is consistent with the typ- 
ical densities where water masers are eventually pumped. The 
upper/lower limits in the number density refer to the B s determi- 
nations for the two limit Ap-F> cases (Table 2). 

If the magnetic field plays a dominant role over the outflow 
dynamics, the magnetic pressure should be similar to the shock 
ram pressure, i. e., 



B 2 J8n : 



: Pov s 



(4) 



where po is the pre-shock density. For the field strength B s deter- 
mined in the previous section and shock velocities of Ristorcelli 
et al. (2005), we find a mean pre-shock density of ~ 5.1 x 10 8 
cm 4 , which is consistent with the pre-shock densities calculated 
from equation 3 and in line with a magnetically controlled out- 
flow evolution. 

Assuming magnetic flux-freezing, the compression of an or- 
dered magnetic field (Bo) amplifies the field strength in propor- 
tion to the gas density as follows 



So 
"o 



(5) 



where n s is the post-shock density. If we consider that the shock 
compresses the well established field geometry in the submm 
maps of Rao et al. (2009), which predict a (scaled) Bo - 9 mG 
at densities no - 4.9 x 10 7 cm 4 , and using the post-shock field 
intensity B s calculated from our data, we find post-shock den- 
sities of ~ 1.3 x 10 9 cm 4 , which is expected for effective wa- 
ter maser pumping. In a more conservative approach, we rescale 
the magnetic field strength of Rao et al. (2009) using B oc rf, 
where k may assume distinct vales depending on the dynamics 
of the core: k = 0.5 for a cloud evolution controlled by the mag- 
netic field (Crutcher et al. 1999), k = 0.47, magnetic evolution as 
modeled by Fiedler & Mouschovias (1993) and k = 0.67, for a 
free-fall cloud collapse. For the three cases, the pre-shock mag- 
netic field in the calculated pre-shock densities ranges between ~ 
24 and 36 mG. Applying again equation 5, we obtain post-shock 
densities ranging between ~ 2 - 5 x 10 9 cm 4 , again consistent 
with typical H2O maser densities, but unable to discern in which 
dynamical regime, if magnetic or turbulent, the core evolves. 

The H2O maser emission of 116293 is strong (> 200 Jy) 
and stable over a few weeks but highly variable over months, 
as shown by surveys performed toward this source (Wilking & 
Claussen 1987; Claussen et al. 1996; Furuya et al. 2003). We are 
then encouraged to carry out VLBI observations in the future. 
With these observations, the blended maser components found 
in this work will be likely resolved both spatially and spectro- 
scopically. We will then be able to apply radiative transfer mod- 
els such as in Vlemmings (2006) to derive the full 3-D magnetic 
field properties as already done in other star-forming regions (e. 
g., Cepheus A HW2, Vlemmings et al. 2006). 

7. Conclusions 

This work reports on the H2O (616 - 523) maser emission ob- 
served with the VLA/EVLA toward the low-mass source IRAS 
16293-2422. This is one of the first estimations of the line-of- 
sight magnetic field strength in a low-mass protostar at densities 
larger than 10 8 cm 4 . From this study, we obtained insights on 
the dynamics of this object at such high densities. Our main con- 
clusions are: 

- We detect strong water maser emission associated with the 
submm source Aa. 

- The maser spectrum has a non-Gaussian profile, indicating 
the presence of at least three maser features covering a veloc- 
ity range of ~ 3.5 km s 4 distributed in a linear configuration. 
Those components are likely excited in zones of compressed 
gas produced by shocks in the outflow activity of 116293. 

- The obtained Stokes V spectrum is consistent with Zeeman 
emission. The mean LOS magnetic field strength is ~ 
113 mG. The POS field topology retrieved from the linear 
polarimetry shows an ordered pattern consistent with larger 
scales field morphologies. 

- The post- and pre-shock densities calculated from our field 
strength estimation are consistent with the typical densities 
expected for H2O maser pumping. 

- The dynamics of the outflow evolution in these sources are 
likely regulated by the magnetic field, since that the magnetic 
pressure in the shocked gas is similar to the pre-shock ram 
pressure. 
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